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Abstract. We present a framework for analyzing luminescence regulation during 
quorum sensing in the bioluminescent bacterium Vibrio harveyi. Using a simplified 
model for signal transduction in the quorum sensing pathway, we identify key 
dimcnsionless parameters that control the system's response. These parameters 
Qi | are estimated using experimental data on luminescence phenotypes for different 

mutant strains. The corresponding model predictions are consistent with results from 
other experiments which did not serve as inputs for determining model parameters. 
Furthermore, the proposed framework leads to novel testable predictions for 
luminescence phenotypes and for responses of the network to different perturbations. 
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1. Introduction 

Bacterial survival critically depends on regulatory networks which integrate multiple 
inputs to implement important cellular decisions. A prominent example is the global 
regulatory network involved in "quorum sensing" , commonly defined as the regulation of 



gene expression in response to cell density During the process of quorum sensing (QS), 
bacteria produce, secrete and detect signaling molecules called autoinducers (Miller 
and Bassler 2001; Waters and Bassler 2005; Bassler and Losick 2006). These signals 
are then processed by the QS pathway to regulate critical bacterial processes such as 
biofilm formation and virulence. The observation that quorum sensing is linked to both 
biofilm formation and virulence factor production suggests that many virulent bacteria 
can be rendered nonpathogenic by the inhibition of their QS pathways (Bjarnsholt and 
Givskov 2007). Quantitative modeling of the QS pathway can thus provide useful inputs 
for treating many common and damaging bacterial infections. 

One of the most studied model organisms for QS based regulation is the 
bioluminescent marine bacterium Vibrio harveyi (Nealson et a!., 1970). Experimental 
studies have led to a detailed characterization of regulatory elements in the pathway 
(Henke and Bassler 2004; Mok et al., 2003; Timmen et al., 2006; Waters and Bassler 
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2006; Tu and Bassler 2007). The network (see Fig. 1) includes multiple autoinducers and 
corresponding sensor proteins which act together to control the phosphorylation of the 
response regulator protein LuxO. The phosphorylated form of LuxO (LuxO-P) activates 
the production of multiple small RNA (sRNA)s which in turn post-transcriptionally 
repress the QS master regulatory protein LuxR. At low cell density, the sRNAs are 
activated and act to effectively repress LuxR expression. In contrast, sRNA production 
is significantly reduced at high cell density, thereby giving rise to increased levels of LuxR 
which leads to the activation of luminescence genes. The corresponding luminescence 
output per cell profile (i.e., colony luminescence/cell output as a function of cell density) 
is frequently used as a reporter to characterize the state of the QS pathway. 

Recent experiments (Henke and Bassler 2004) have analyzed the effects of 
mutagenesis of different pathway components on the corresponding luminescence profile 
in V. harveyi. It was observed that there are distinct luminescence profiles as the 
network is perturbed corresponding to different pathway mutants. The changes in the 
luminescence profile were used to infer pathway characteristics such as relative kinase 
strengths for the different sensors. Given the complexity of the network which involves 
integration of multiple inputs, it would be desirable to develop a quantitative framework 
for inferring pathway characteristics based on network perturbations. The corresponding 
quantitative model can then be used to make testable predictions for future experiments 
as well as to further analyze existing experimental data. The aim of this work is to 
develop such a minimal model for the QS pathway in V. harveyi. 

The starting point of our analysis is the observation that luminescence/cell output 
is controlled by the degree of phosphorylation of the response regulator LuxO. We thus 
develop a simplified model which connects external autoinducer concentrations to the 
degree of phosphorylation of LuxO for the wild type (WT) strain and for different 
mutants. Our analysis identifies key dimensionless parameters which control the system 
response and which can be determined using the experimental results for luminescence 
phenotypes. Determination of the effective parameters, in turn, leads to predictions 
for the systems response to a broader range of perturbations, i.e., perturbations distinct 
from those used to infer the effective parameters. The corresponding analysis sheds light 
on previously obtained experimental results and also gives rise to testable predictions 
for future experiments. 

The rest of the paper is organized as follows. In Section 2 we give an overview of 
the QS network in V. harveyi. We then develop a minimal model of the QS pathway and 
define key dimensionless parameters which control the network response characteristics. 
In Section 3, we connect our model to experimental data on different luminescence 
curves and thereby determine model parameters. In Section 4, we discuss experimentally 
testable predictions based on the model and conclude with a summary. 
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Figure 1. (color online) Schematic representation of quorum sensing network in Vibrio 
harveyi at high and low cell densities. The dotted rectangle is the input module which 
controls phosphorylation of LuxO in response to external autoinducer concentrations. 
Solid line, active path; Dotted line, inactive path; IM, inner membrane; OM, outer 
membrane. 



2. Overview and Model 



The QS network in V. harveyi is shown in Fig. 1. The key upstream components of 
the pathway are the three sensors, LuxN, LuxPQ and CqsSy^ and the corresponding 
autoinducer synthases, LuxM, LuxS, and CqsAy^ which are responsible for producing 
the three autoinducers: H-AI1, AI-2, and CAI-1, respectively. The binding of a 
single autoinducer to a sensor is highly specific, i.e., HAI-1 binds only to LuxN, AI- 
2 binds to LuxPQ only, and CAI-1 binds specifically to CqsSy/j (see Fig. 1). The 
overall network is conveniently described in terms of functional modules. The first 
(input) module includes interactions between autoinducers ([ALJ (i = 1,2,3)) and the 
corresponding sensor proteins which, through a phosphorelay mechanism, determine the 
overall phosphorylation state of a <r 54 -dependent response regulator LuxO. 

The second module focuses on the regulated production of sRNAs (dependent 
on the phosphorylation state of LuxO) and the interaction between the sRNAs and 
the master regulator protein, LuxR. The interactions between small RNAs and their 



A model for signal transduction during quorum sensing in Vibrio harveyi 4 




Cell Density 

Figure 2. (color online) Schematic representation of typical luminescence curves 
from experiment. The green curve represents the response from a wild type (WT) 
colony. The turnaround point in the curve corresponds to cell density necessary for 
the activation of the genes responsible for luminescence output per cell. The red curve 
represents the luminescence/cell curve for a mutant strain that is able to achieve the 
same activation at a lower cell density. 

regulated targets have been modeled in several recent studies which shed light on how 
target protein expression is controlled by small RNA-mediated regulation (Lenz et al., 
2004; Levine et al., 2007; Levine and Hwa 2008; Mehta et al., 2008; Mitarai et al., 
2007). In V. harveyi, LuxR serves as the target protein whose expression is controlled 
by the small RNAs in combination with the RNA-binding protein Hfq. The resulting 
concentration of LuxR determines the level of activation or repression of a multitude 
of genes including the genes involved in bioluminescence (Waters and Bassler 2006). 
The corresponding change in the luminescence/cell output determines the luminescence 
profile which is frequently used to infer network characteristics such as relative rates of 
kinase/phosphotase activities by the sensor proteins (Henke and Bassler 2004). 

A schematic representation of typical luminescence/cell curves is shown in Fig. 2. 
Since the starting point is obtained by the dilution of cells in the high density limit, the 
luminescence output per cell is maximal at the initial time points. The luminescence 
output per cell then declines steadily with increasing cell density, since luminescence 
genes are no longer activated in the cells. At a specific cell density, the luminescence 
curve starts to rise again signalling the start of de novo luminescence gene activation 
by cells in the growing colony. The cell density necessary for activation can vary from 
the WT and mutant strains resulting in different luminescence phenotypes (see Fig. 2). 

Current data indicates that increasing cell density leads to increasing dephosphory- 
lation of LuxO leading to lower production rates for the sRNAs. Correspondingly, the 
turnaround point in the luminescence curves corresponds to unphosphorylated LuxO 
reaching a critical level above which sRNA production is not effective at repressing 
LuxR levels below the threshold for observable luminescence activation in the popu- 
lation of cells. Thus, understanding how external signals (i.e., Al concentrations as a 
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function of cell density) are translated into the degree of LuxO phosphorylation (i.e., 
the input module) is critical for analyzing luminescence profiles. Furthermore, pathway 
mutants which function upstream of LuxO are not known to have any direct effects on 
sRNA production or LuxR levels, apart from the indirect effects mediated by LuxO. 
Therefore we expect that the critical level of LuxO phosphorylation corresponding to 
the turnaround in the luminescence profile is the same for all mutants. The observa- 
tion that the luminescence profiles are different for different pathway mutants indicates 
different functional relations between external AI concentrations and LuxO phosphory- 
lation levels for the different mutants. In the following, we derive a simple model which 
connects cell density to LuxO phosphorylation and uses information from luminescence 
profiles of different mutants to infer system parameters. 

The sensor proteins in the QS pathway can be modeled as two state systems 
(Neiditch et al., 2006; Swem et al., 2008). We consider a further simplification which 
takes the sensors to be existing either in the kinase mode, S^, or in the phosphatase 
mode, S P i (where i — 1,2,3 corresponds to the distinct sensor proteins in V. harveyi). 
In the kinase mode, the sensors can autophosphorylate and then transfer the phosphate 
group to the downstream protein LuxU, whereas in the phosphatase mode the phosphate 
flow is reversed. Experiments indicate that at low cell density (corresponding to low 
autoinducer concentrations) the sensors are primarily in the kinase mode, whereas at 
high cell density (corresponding to high autoinducer concentrations), the sensors are 
primarily in the phosphatase mode. Correspondingly, we consider a simplified model 
wherein the free sensor corresponds to the kinase mode, whereas binding of autoinducer 
results in a transition to the phosphatase mode. 

At a given cell density, the external autoinducer concentrations will be proportional 
to the colony forming units N. Since the time scale for changes in N (i.e., the 
doubling time) is large compared to the time scales for binding/unbinding of ligands 
and subsequent phophorylation/dephosphorylation, the corresponding reactions can be 
considered in steady state for a given N. Furthermore, since the typical number of sensor 
proteins of each type is large, the concentration of sensors of type % is well approximated 
by the mean value [Si] = q[£o] (where [So] is some reference concentration). At a given 
cell density, external AI concentrations determine the fraction of the receptors which 
exist in either the kinase or phosphatase mode. For the simplest case of autoinducers 
binding to their cognate sensors, we have the kinetic scheme: 



from which the mean steady state concentrations of the sensors in either the kinase or 
phosphatase mode can be obtained. More generally, to account for cooperative effects 
in binding, we take the kinase/phosphatase fractions to be: 



ki + All =F=^ Dpi, 
k—i 



(1) 



[Ski] = (1 - 9i)ci[So] and [S pi ] = giCi[S ], 



(2) 



where 



[Ski] + [Spi] = Ci[S ], gi = </(l + <), a, 



[AIA/Ki. 



(3) 
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and hZi — fc—i I 

Equation (j2J), with Hill coefficient n — 1, corresponds to the steady state fractions for 
equation (JTJ), higher n values correspond to sharper switching from kinase to phosphatase 
mode which mimics cooperative effects in binding. Finally, since the concentration 
of the i-th autoinducer, [AZj], is proportional to the colony forming units (CFU), N, 
i.e. [Ali] = ViN] we renormalize the binding constant Ki to define the scaled effective 
parameter Ki = Ki/vi. 

Typically in bacterial signal transduction, the sensor proteins in the ki- 
nase/phosphatase modes serve as enzymes which transfer the phosphate group to/from 
a response regulator protein or a phosphorelay protein (Appleby 1996; Hoch 2000; Stock 
et al., 2000; Laub and Goulian 2007). In V. harveyi, this step involves phosphotransfer 
to the phosphorelay protein LuxU (U). Phosphorylated LuxU (Up) can then transfer 
the phosphate group to the response regulator LuxO (O); similarly, unphosphorylated 
LuxU serves as a receiver for removing the phosphate group from phosphorylated LuxO 
(Op) . We represent these processes by the following equations: 

S ki + U h -^S ki + Up, (4a) 

S pi + U P ^ S pi + U, (46) 

Up + O %U + P . (4c) 

k b 

For the above kinetic equations, it is convenient to define key dimensionless 
parameters of the model as follows 

an = Cik ki /k kr} j3i = (kf/k b )(k ki /k P i). (5) 

The parameter a r i is a measure of the relative kinase strength of i-th sensor with respect 
to the r-th sensor (scaled by the mean concentrations of the two sensors), e.g., a±2 is 
the relative kinase strength of sensor 2 with respect to sensor 1. Another set of key 
parameters is the ratio of the scaled kinase to phosphatase rates, of the i-th sensor. 
Using these dimensionless parameters, we then solve the rate equations (4a-4c) at steady 
state to derive the following expression for the fraction of unphosphorylated LuxO at 
steady state, /luxO = [0]/[0]o (with [0]q being the total LuxO concentration) 

f L Q = E»«»fa/A) , Q s 

3. Connection to experimental data 

We now connect the model for LuxO phosphorylation developed in the previous section 
to experimental luminescence curves. Recall that the typical luminescence profile 
shows a well defined switching point which signals observable de novo production of 
luminescence by the population of cells. As argued earlier, this corresponds to a critical 
value for the concentration of unphosphorylated LuxO. Let us denote this critical 
fraction of unphosphorylated LuxO by f c and the corresponding value of the colony 
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forming units by N c . At /luxO — / c ? f° r the WT luminescence curve we have the 
following relation ([6]): 



where the factors gi are evaluated at N = N c . Since N c is known from experiments 
corresponding to the WT luminescence curve, the above equation can be regarded as a 
constraint on the dimensionless parameters. 

We now consider the corresponding equations for luminescence phenotypes 
of the mutant strains. Current knowledge of the QS network in V. harveyi 
indicates that pathway proteins functioning upstream of LuxO primarily control LuxO 
phosphorylation levels and have no direct interactions with the qrr sRNAs or the master 
regulator LuxR. This suggests that for each mutant the degree of LuxO phosphorylation 
needed to activate luminescence is the same (i.e., f c is the same) since upstream 
proteins affect LuxR only via LuxO-P levels. The observation that the luminescence 
profiles are distinct for different pathway mutants is a consequence of the altered 
functional relationship between LuxO phosphorylation levels and external autoinducer 
concentrations for the mutants. Given the defined roles of the pathway proteins, these 
altered functional relationships can readily be derived within our model for all the 
mutants. For example, equation © for the single sensor mutant cqsSvh (i-e. the strain 
with a deletion for the gene cqsSvh) takes the form: 



Note that the quantity (1 — f c ) / f c can be absorbed into the scaled kinase to phosphatase 
ratios (3i and f3 2 - This is equivalent to setting f c = 1/2 in the above equation, and 
since f c is the same for all pathway mutants, a similar rescaling can be done for the 
functional relationships for all the mutants. The corresponding equations are presented 
in the Appendix. In the following, we show how these equations can be used along with 
WT and mutant luminescence phenotypes to determine effective system parameters and 
to make testable predictions. 

From the work of Henke and Bassler (2004), the critical threshold in colony forming 
units (N c ) can be estimated for a range of pathway mutants. The different mutant 
strains studied were i) luxN, ii) luxQ, iii) cqsSyh, i y ) luxN luxQ, v) luxN cqsSyh, and 
vi) luxQ cqsSvh- To connect the sensors of V. harveyi with our model, we designate 
sensors LuxN, LuxQ, and CqsSy^ as 1, 2, and 3, respectively. The ordering of the 
CFU/volume for the different strains at their critical threshold shows the following 
hierarchy (Henke and Bassler 2004): 



where N£ 2 is the number of colony forming units for mutant strain luxN luxQ at 
which /luxO = f c and so on. Although the values iV|, N 23 and iVf, iVf 3 appear to be 
indistinguishable based on available experimental data, based on the model developed 




(7) 




N c l2 « N c 2 ~ N C 2Z < N c wt < iV 3 c < N c x 



(8) 
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we expect a small difference in the threshold values. For example, the difference between 
the luxN strain and luxN cqsSyh strain is that CqsSy^ is active as phosphatase in the 
luxN mutant (close to the switching threshold). This implies that the switching in the 
luminescence phenotype should occur at a lower N c value for the luxN cqsSyh strain i.e., 
iVf < N£ 3 . Since CqsSy^ has weak effect on the luminescence phenotype, the switching 
values are indistinguishable experimentally. However to develop a consistent model, we 
have to impose a small difference between the switching values based on the constraint 
iVf < Ni 3 (and similarly for N 2 and N 23 ). 

Based on the above reasoning, we initially considered a ~10% difference between 
N 2 , N 23 and iVf, N£ 3 . Accordingly, the values for critical thresholds (switching values, 
in the units of CFU/volume) used as initial inputs were 

N{ 2 ~ 1 x 10 5 , iV 2 c ~ 14 x 10 5 , iV 2 C 3 ~ 15 x 10 5 , N c wt ~ 40 x 10 5 , 
iV 3 c ~ 70 x 10 5 , ~ 110 x 10 5 , N£ ~ 100 x 10 5 . 

From the discussion of the previous section, we have seen that the input module 
provides us eight key parameters: two relative kinase strengths (a i2 and a^), three 
scaled kinase to phosphatase ratios (fli, fa, and $3) and three effective binding constants 
(ki, k 2 , and k 3 ). Given that we have experimental data for threshold cell densities for 
seven strains, this indicates that if one of the parameters is fixed, the other parameters 
can potentially be determined by solving the corresponding threshold equations (see 
Appendix). Since previous work indicated that the effect of CqsSyh on luminescence 
phenotypes is minimal, we initially fixed the parameter a±3 (the relative kinase strength 
of sensor 3 (CqsSy/J with respect to sensor 1 (LuxN)) to 0.0010. We then proceeded to 
determine the effective model parameters by solving the threshold equations using the 
above experimental inputs for switching cell densities. We also checked the stability 
of the solutions to the above equations based on small perturbations to the input 
parameters (data not shown). We found that the solutions are stable with respect 
to perturbations that maintain the initial ~10% difference between N 2 , N 23 and N£, 
N± 3 . However the solutions are sensitive to changes in the parameters controlling the 
small differences in N c values. Since experiments cannot guide us in determining the 
precise value of these differences, the values of N 2 and N° do not serve as useful inputs 
in determining model parameters. Thus additional experimental data is needed to 
determine model parameters as outlined below. 

The luminescence data at high cell densities (hcd) for different sensor mutants 
from the work of Henke and Bassler (2004) (see Fig. 4A) provides an indirect means of 
estimating model parameters. The basic experimental observations can be summarized 
as follows: while the WT strain shows a bright phenotype at hcd, the luxS strain 
has a dim phenotype and the luxM strain has low levels of luminescence and is 
classified as being dark. Furthermore the cqsSyh strain has a luminescence output 
that is intermediate between WT and luxS and the cqsA Vh luxS double mutant is 
dark and produces significantly less luminescence than a luxM strain. Given our 



% This assumption will be relaxed in the subsequent analysis as described below 
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colony forming units (CFU)/volume 

Figure 3. (color online) Profile of /l U xO as a function of colony forming units 
(CFU)/volume for wild type (WT) and different sensor mutant phenotypes. The cell 
density at which /luxO = f c corresponds to the turnaround point in the experimental 
luminescence curves. 



definitions of model parameters, /luxO = 1/2 corresponds to value at which observable 
luminescence/cell is produced. Higher values of /luxO will correspond to brighter 
luminescence phenotypes, whereas a dark luminescence phenotype implies /luxO < 1/2. 
Thus we expect that, at hed, we have /luxO f° r luxS mutants to be around 0.5 (given the 
dim luminescence phenotype) and /luxO for the cqsSyh strain to be significantly greater 
than the corresponding value for the luxS strain but significantly lower than 1 (the value 
for the WT strain). Based on these constraints, we set the /luxO values for 3 synthase 
mutant strains at hed as follows: /luxO = 3/4, /[uxo = an d /luxo"^ = 1/4. In 
combination with the expression derived for /luxO (equation (J6])), these equations can 
be used, along with luminescence switching cell density equations, to determine model 
parameters (see Appendix). 

First, considering equation (CO) for the double sensor mutants, we have the relation 
between the three /3-s and three R-s, 

AT 2 C 3 = R^\' n , N c 13 = K 2 /3 2 1/n , N c 12 = R 3 (3l /n . (9) 
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Also from equation ([7]), we have the expressions for the wild type and one single sensor 
mutant (cqsSvh) with five unknown parameters: three kinase to phosphatase ratios 
(Pi, Pi and Ps) and two relative kinase strength (au and a^) (Note that we are now 
considering a 13 to be variable). Using /luxO — f c — 1/2, the switching values for WT 
and cqsSvh, (N^ t and N£) and the /luxO for three synthase mutants at hcd we solve 
the five equations to determine the five unknown parameters. The corresponding values 
for the key parameters of the model are: cti 2 ~ 0.14, a i3 ~ 0.19, Pi ~ 8.99, P 2 ~ 0.29 
and Ps ~ 7.14, for the Hill coefficient n = 4. We note that there are two sets of 
solutions obtained using the above approach, however only one of these corresponds to 
the experimentally observed hierarchy of switching cell densities flH]). Furthermore no 
solutions were obtained for n ^ 2. For n — 3, the equations can be solved and yield 
parameters that are close to the those inferred for n = 4. However the n = 4 results are 
more consistent with the experimental observation that the switching cell densities are 
experimentally indistinguishable for N£ and N£ 3 (similarly for iV| and iVJg). The high 
value of n = 4 suggests that there might be cooperative effects in the switch from the 
kinase to phosphatase mode for the sensors. Now using these values for the effective 
parameters, we calculate the values of /l ux o as a function of CFU/volume (see Fig. 3) 
for the WT and different sensor mutant phenotypes of V. harveyi. Since the effective 
parameters are determined, we can now use our model to generate similar curves and 
make predictions for mutants that have not yet been studied experimentally. We have 
checked the stability of the obtained solutions with respect to small changes in the input 
values (see Appendix). We have also considered larger changes in the input /luxO values 
consistent with the constraints noted earlier. While the precise values of the effective 
model parameters do change as the inputs are varied, there are several robust predictions 
that can be made. These are discussed further in the concluding section. 

4. Conclusion and outlook 

The preceding analysis helps determine the parameters in our minimal model. While 
these parameters cannot directly be compared to experiments, they can lead to several 
predictions which are testable experimentally. In the following, we outline some of the 
key predictions based on our analysis. 

1) The parameter Pi is a measure of the relative kinase to phosphatase rates for 
the i-th sensor. Based on the values determined, the following ordering is predicted for 
the relative kinase to phosphatase rates of the three sensors LuxN > CqsSy^ >LuxQ. 
LuxN is predicted to be the strongest kinase which is consistent with results from 
previous experiments showing that LuxN has a greater effect on LuxO phosphorylation 
than LuxQ (Freeman et al., 2000). Furthermore, it is interesting to note that recent 
experiments have demonstrated high kinase to phosphatase rates for the sensor LuxN 
(Timmen et al., 2006). While the corresponding value estimated by our model (Pi ~ 9) 
cannot directly be compared to experiments since it involves additional parameters, the 
ratio Pi/ Pj (i ^ j) should correspond to experimental estimation of the ratio of kinase to 
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Table 1. Predictions for luminescence output per cell of different synthase mutants 



and mixed 


sensor- synthase mutants. 


Phenotype 


Mutant 


dark 


luxM, luxM luxS, luxS cqsAyh, luxM cqsAyh, 




luxN luxS, luxQ luxM, luxQ cqsAyh, cqsSyh luxM 


dim 


luxS, cqsSvh luxS 


bright 


cqsAyh, luxN cqsAyh 



phosphatase rates of two sensors. From our model we consistently find that j3 2 / Pi <C 1 
and /3 2 //3 3 <C 1 indicating the the effective kinase to phosphatase activity ratio for LuxQ 
is much lower than the other two sensors. Note that this prediction differs significantly 
from the previous characterization (Henke and Bassler 2004) that kinase to phosphatase 
activity ratio for LuxQ is greater than that of CqsSyh- It would thus be of interest to 
carry out experiments to measure relative kinase to phosphatase rates for the sensors 
LuxQ and CqsSyh to see if the predictions are borne out. 

2) Experiments with mutant strains (besides those used as inputs to our model) 
indicate that at high cell densities, the luminescence phenotypes can be broadly 
categorized into 3 types: dark, dim and bright. Since f c = 1/2 is the threshold 
for luminescence activation in our model, we take these categories to correspond 
to the following: dark (0 ^ AuxO < 0.4), dim (0.4 < /luxO < 0.6) and bright 
(0.6 < /luxO ^ 1-0). Using these criteria, we can now predict the luminescence 
phenotypes at high cell density for other pathway mutants (i.e. those not included in the 
experimental inputs used to determine model parameters). The corresponding results 
are listed in Table [lj We note that all mutant strains with LuxM deleted (luxM) are 
dark. This is consistent with previous experimental results (Freeman and Bassler 1999). 
Other interesting predictions are 

i) While cqsA Vh luxN is bright (comparable to cqsA Vh ) at hed, the strain cqsA Vh luxQ 
is predicted to be dark; 

ii) luxS is brighter than luxM at hed , however cqsA Vh luxS is predicted to be darker 
than cqsAyh luxM (note that this is consistent with the observations in Henke and 
Bassler (2004)). 

It should be noted that the results presented in Fig. 3 are just for sensor mutants 
whereas Table [1] is for synthase mutants and mixed sensor-synthase mutants. For the 
different mutants given in Table [U, the maximal value of the /luxO curve differs from 
1 and stays within the defined range (according to the broad categories discussed in 
the paper) even at the hed in contrast to the behavior shown in Fig. 3 for the sensor 
mutants. 

3) To figure out the values of the effective parameters of the model, we have used 
the switching value {N c ) of WT, cqsSyh and double sensor mutants from the experiment 
(Henke and Bassler 2004). With these derived values of the effective parameters, we can 
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now predict the switching values of the other two bright sensor mutant strains (luxN 
and luxQ) at hcd (in the units of CFU/volume), 

N{ ~ 100 x 10 5 , iV 2 c ~ 14 x 10 5 . 

It is interesting to note that the above switching values are in good agreement 
with the observation that N{ is experimentally indistinguishable from N° 3 and N% 
is experimentally indistinguishable from N% 3 (see Fig. 3). In addition, the effective 
parameter set predicts the switching values (N c , in units of CFU/volume) for the two 
bright mutant strains cqsAyh and luxN cqsAyh mentioned in Table [Has ~ 130 x 10 5 
and ~ 156 x 10 5 , respectively. 

4) Recent experiments have probed the response of the QS pathway to externally 
controlled autoinducer concentrations (Mok et al., 2003). In these experiments, the 
autoinducer production is switched off by deleting the corresponding synthases and 
then autoinducers are added back exogenously in controlled amounts. In our model 
this behavior can be mimicked by controlling the quantity in equation ([3]). For 
each synthase mutation the autoinducer production is switched off so that gi — as 
Ali — (i — 1,2,3). As autoinducers are added to the network from outside, the 
quantity g^ grows and tends to one as Ali °o. For this setup, our analysis indicates a 
situation wherein the sensor CqsSy^ plays an important role in regulating the response 
which is contrary to what is normally assumed. Consider the situation for which all 
the autoinducer synthases have been deleted and subsequently saturating amounts of 
Ali are added. In this case, we predict a significant difference between the luminescence 
output per cell for the two cases corresponding to i) low external AI 3 concentrations and 
(ii) high external AI3 concentrations. The difference between these two cases is that the 
sensor CqsSy^ is primarily in kinase mode for case (i) and in phosphatase mode for case 

(ii) . Our analysis thus suggests a testable prediction for an experimentally realizable 
situation wherein signaling through CqsSy^ significantly changes the output from the 
QS pathway. 

5) Finally, we examine predictions from our model for the expression of genes 
that are also controlled by /luxO through LuxR but are not directly related to 
luminescence/cell. Waters and Bassler (2006) studied several genes regulated by LuxR 
and classified them into different categories based on the activation/repression induced 
by the presence of high concentrations of either Ali or AI2 or both. We will focus on 
the category of genes (labeled "class 3" genes) which are defined as genes that show an 
equally notable change in expression when either Ali and/or AI2 are present in high 
concentrations. Within our model, we can calculate the the values of /luxO for the 
3 cases : (i) High concentration of AI X only, (ii) high concentration of AI 2 only and 

(iii) high concentration of both Ali and AI 2 . Out of these the lowest value of /luxO 
corresponds to case (ii) i.e., high concentration of AI 2 only. Since class 3 genes are fully 
activated/repressed when high concentrations of AI2 only are present, it follows that 
the f c for all genes in this category must be lesser than the value of /luxO when only 
AI 2 levels are high (/luxO = 0.33). (Note that we have assumed that AI 3 levels are at 
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high concentrations in the above experiments since they are at high cell densities). This 
observation indicates that an upper bound for activation/repression of class 3 genes 
corresponds to f c = 0.33. Using this, the following testable predictions can be made 

• The synthase mutant luxM can fully activate/repress class 3 genes at high cell 
density. Note that luminescence genes, in contrast, are not activated at high cell 
density in a luxM mutant. 

• Similarly, the sensor-synthase mutants luxM cqsSyh and luxQ cqsA Vh cannot 
activate luminescence genes at high cell density whereas they are predicted to fully 
activate/repress all class 3 genes at high cell density . 

The minimal model presented in this work can be generalized further as more 
experimental data becomes available. An important generalization would be to relax 
some of the assumptions made by considering a two-state model (Swem et al., 2008) 
which incorporates non-zero phosphatase activity in the on (free) state and nonzero 
kinase activity in the off (bound) state. We note that this will add several additional 
parameters to our current model. With additional experimental data, the generalized 
model could be used to estimate the expanded set of effective parameters. While 
the effective parameters so determined are likely to be different from the values 
determined using the minimal model, the framework connecting the model parameters 
to experimental data will essentially be the same. 

In summary, we have proposed a minimal model to study the quorum sensing 
network in V. harveyi. Using experimental data for luminescence phenotypes of WT and 
different mutant strains, we provide a framework to estimate the effective dimensionless 
parameters of the model. Correspondingly, the model can be used to predict the 
luminescence phenotypes of other pathway mutants which have not been experimentally 
studied to date. The proposed framework captures the key features of the signal 
transduction in V. harveyi and can contribute to guiding and interpreting experimental 
efforts analyzing the QS pathway in the Vibrios. 
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Appendix 

For the relative kinase strength (a r i = Cikkt/kk r for % = 1,2,3) of the sensors, we 
generally use the kinase strength of LuxN, i.e., kki(r = 1), as the reference kinase. Now 
using equation (JTj) we explicitly write the functional relation for WT strain evaluated 



at N = N° wt for / Lux0 = f c : 



(1 - g x ) + a 12 (l - g 2 ) + «i3(l - 9z) 



( 



l-f c \ \9i 



f c ) m 




(A.l) 
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Similarly for luxN mutants we use kinase strength of LuxQ, i.e., k^ 2 [r = 2), as the 
reference kinase whereas for luxQ and cqsSvh we use kinase strength of LuxN as the 
reference kinase as in WT. Thus the functional relations for the single sensor mutants 
luxN, luxQ and cqsSyh evaluated at N£, N 2 and iV|, respectively, are: 
For luxN (r=2): 

1 /' 



;i -92) + — (l-ss) 

"12 



For luxQ (r=l) 



1 - gi) + ai 3 (l - g 3 ) 



For cqsSyh (r=l) 
fl- 



9i) + «i 2 (l - g 2 ) 



f 

1-f 
f 

1-f 



^2 «13^3 
,/?2 "12 As 



01 . #3 



01 . #2 



(A.2) 



(A.3) 



(A.4) 



For double sensor mutants, the value of the relative kinase strengths becomes 1 as there is 
only 1 sensor. Hence the functional relations for the double sensor mutants luxN luxQ, 
luxQ cqsSyh and luxN cqsSyh evaluated at iVf 2 , N 23 and N£ 3 , respectively, are: 
For luxN luxQ (r=3): 



(i-o») = 

For luxQ cqsSyh ( r= l) : 
(l-<7i) = 
For luxN cqsSyh ( r =2) 

(1 " 92) = 



1-f 
f 



1-f 

f 



1-f 

f 



93 



£l 
A' 



02 



(A.5) 



(A.6) 



(A.7) 



Now, using equation (Q and the expression for the functional relation for the WT strain 
given above (A.l) we have the following equation evaluated at N — N° t for f c = 1/2, 

1 + {N^/RtY 



1 - 



+ «13 



+ «12 



1 + (N< t /R 2 ) 



1 + (N^/R.r 
1 + {N'jRzY 



1 - 



N: 



wt 



N: 



12 



0. 



Similarly, for luxN, luxQ and cqsSyh strains we have the following set of equations 
evaluated at N°, N 2 and iV|, respectively, for f c = 1/2, 



1 - 



^23 

Nl 
N c 2 z 



+ 01 



Olg 1 + (iVf /K 2 )" 

aia 1 + (iVf/« 3 ) n 

i + (iv 2 c M) w 

1 + (A^)" 
1 + (iV|/Ki) W 

i + (iv 3 7K 2 )" 



+ 



+ «13 



12 



1 - 



1 - 



12 
^2 

A 3 C 



0. 



13 



(A.9) 
(A.10) 
(A.ll) 
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To find the unknown parameters of the system of equations (ai 2 , «i3, Pi, P2, and P3), 
we use equations (1A.8|) and flA.llj) evaluated at N = N° t and A" = N£, respectively, 
along with the following three equations all evaluated at N = N Urge : 



rluxM = (Xl2(92/ P2) + 013(^3/^3) 

/Lux0 l + a 12 (l-g 2 ) + a 13 (l-g 3 ) + ai 2 (g 2 /p 2 ) + a 13 (g 3 /p 3 y { ' ' 

rcqsA = (9l/ Pi) + <Xl2(g2/ fo) / a 1 o\ 

/Lux0 ai 3 + (l-gi) + a 12 (l-g 2 ) + (g 1 /p 1 ) + ai2(g2/P2y 1 ' ' 

AuxS cqsA _ (9l/ Pi) ,^ A , 

/Lux ° a 12 + a 13 + (1 - g x ) + (g x /p x ) ' 1 ' ] 

Equations (1A. 12HA. 14|) give the /luxO values for the three mutants luxM, cqsA V h, and 
luxS cqsAvh once the system has reached steady-state for N = N Urge such that /luxO 
has saturated. Equations (1A.8I) and (lA.lltiA.lll) are then numerically solved using 
Mathematica (Wolfram Research, Inc., Version 6, 2008) which yielded two solutions 
subject to the constraint that all the parameters must be real and positive. We keep the 
solution that best agrees with experimental data (see main text). When solving these 
equations, we used f*£% = 0.33, = 0.75, f l ^ Q cqsA = 0.25 and n = 4. 

We next analyzed the changes to the solutions based on small perturbations to 
the input parameters. Each perturbation for the input values is drawn from a random 
Gaussian distribution whose mean is the base value and variance is the base value xa, 
where a is chosen such that 68% (98%) of the perturbed values lie within 2% (5%) of 
the base value. For example, to generate a list of perturbed N{ 2 values, we set the 
mean of the Gaussian distribution to be N± 2 and the variance to be N± 2 x o~, etc. Using 
this scheme, we generated 100 random data points for the input values (the switching 
values) and numerically solve equations (1A.8I) and flA.lltiA.T4|) with n = 4 to generate 
the effective parameters. Note, fj^/o, /luxO> /luxO^" 4 are a ^ so perturbed in the same 
fashion. 

The resultant data of the sensitivity analysis are shown in Figs. A1-A2. The nature 
of the data shown in Figs. A1-A2 suggests that the parameter set obtained using the 
experimental switching values from Henke and Bassler (2004) is robust against small 
perturbations. 



Glossary 

Quorum sensing. Process leading to regulation of gene expression in response to cell 
density. 

Autoinducers. Small signaling molecules produced by bacteria which bind to specific 
receptors and induce the quorum sensing response. 
Kinase. Enzyme acting as phosphate donor. 
Phosphatase. Enzyme acting as phosphate acceptor. 

Bioluminescence. Production of light by living organism as a result of internal chemical 
reactions. 

Vibrio harveyi. Gram-negative and bioluminescent marine bacterium. 
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Figure Al. (color online) Results of sensitivity analysis for the input base values 
The blue line represents the unperturbed data and the red dashed line is the mean of 
the 100 perturbed data points represented by scattered red points. 
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Figure A2. (color online) Results of sensitivity analysis for the effective parameters. 
The blue line, red line and red scattered points have the same meaning as in Fig. Al. 
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